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The radius and mass of the subgiant star (3 Hyi from interferometry 
and asteroseismology 
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ABSTRACT 

We have used the Sydney University Stellar Interferometer (SUSI) to measure the angular 
diameter of f3 Hydri. This star is a nearby G2 subgiant whose mean density was recently 
measured with high precision using asteroseismology. We determine the radius and effective 
temperature of the star to be 1.814 ± 0.017 R© (0.9%) and 5872 ± 44K (0.7%) respectively. 
By combining this value with the mean density, as estimated from asteroseismology, we make 
a direct estimate of the stellar mass. We find a value of 1.07 ± 0.03 M© (2.8%), which agrees 
with published estimates based on fitting in the H-R diagram, but has much higher precision. 
These results place valuable constraints on theoretical models of (3 Hyi and its oscillation 
frequencies. 
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1 INTRODUCTION 



The combination of interferometry and asteroseism ology is a pow- 
erful way of constraining t he parameters of star s jKervella et al.l 
20031 : IPiipers et al.1 120031; iKervella et al.l |2004 iTheyenin et al.l 



ature and surface gravity are also constrained using our radius and 
complementary data from the literature. 



200l |yan Belle, Ciardi, & Boden 2007: ICreevev et alj|2007l) . Ir^ 2 OBSERVATIONS AND DATA REDUCTION 



terferometry provides an angular diameter which, combined with 
the parallax, gives a direct measurement of the stellar radius. As- 
teroseismology, on the other hand, uses the oscillation frequen- 
cies of the star to infer details about its internal structure (e.g., 
iBrown & Gillilandl 19941) . In particular, Ai/, the so-called large fre- 
quency separation between consecutive radial overtones, gives a 
good estimate of the mean stellar density. The combination of the 
radius and Av can therefore provide a direct measurement of the 
mass of a star. 

The star (3 Hyi (HR 98, HD 2151, HIP 2021) is a southern 
evolved subgiant (spectral type G2 IV) with a mass slightly greater 
than the Sun and an age of about 6.5-7.0 Gy i lDravins et all 19981 : 
iFernandes & Monteiroll20o"3h . This star is an excellent target for 
asteroseismology, being the closest G-type subgiant. Its distance, 
thanks to Hipparcos, is known to 0.36% accuracy. It is bright (V = 
2.8) and has a low vsini, which allows extremely good Doppler 
pr ecision. Solar-l i ke os cillations in [3 Hyi were recently measured 
by Bedding et al] d2007l) . who used the large frequency separation 
to infer the mean stellar density to an accuracy of 0.6%. 

In this paper we report the first measurement of the angular di- 
ameter of (3 Hyi. We use the Hipparcos parallax to infer the radius, 
which we combine with the mean density, determined from astero- 
seismology, to estimate the mass of the star. The effective temper- 



The S ydney University Stellar Interferometer (SUSI, iDavis et all 
1 19991) was used to measure the squared visibility (i.e. normalised 
squared modulus of the complex visibility) or V 2 on a total of 
7 nights. The red- table beam-combination system was employed 
with a filter of centre wavelength and full-width half-maximum 
700 ran and 80 ran respectively. This system is to be described in 
gr eater detail by Davis et al. (in preparation) and an outline is given 
bv lTuthilletai1 ( l2004l) . 

The interference produced by a pupil-plane beam-combiner 
was modulated by repeatedly scanning the optical delay about the 
white light fringe position. Two avalanche photo-diodes detect the 
two outputs of the combining beamsplitter. An observation con- 
sisted of a set of 1000 scans, each traversing 140 /im in optical 
delay, sampled in 1024 steps of 0.2 ms duration. 

The post-processing software mitigates the effect of scintilla- 
tion in each observation by differencing the two recorded signals. 
The squared visibility was estimated after bias subtraction and a 
nonlinear f unction was ap plied to partially correct for residual see- 
ing effects ((Ireland 2006). 

Target observations were interleaved with calibration mea- 
surements of nearby stars chosen to provide a known or essentially 
unresolved signal. The angular diameters (and associated error) 
of these calibrators stars were estimated from an intrinsic colour 
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Table 1. Adopted parameters of calibrators used during observations. 



HR 


Name 


Spectral 


V 


UD Diameter 


separation 






Type 




(mas) 


from /3 Hyi 


0705 


<5Hyi 


A3V 


4.08 


0.50 ±0.03 


11?85 


0806 


eHyi 


B9V 


4.10 


0.39 ±0.03 


13? 05 


7590 


e Pav 


A0V 


3.96 


0.46 ±0.05 


16?53 



interpolation (and spread in data) of mea surements made with 
the N arrabri Stellar Intensity Interferometer jrlanburv Brown et al.l 
1 19741) . including corrections for the effects of limb-darkening. Us- 
ing the adopted stellar parameters of the calibrator stars given in 
Table Q] the system response or transfer function was quantified. 
The weighted mean of the bracketing transfer functions was then 
used to scale the observed target squared visibility appropriately 
to produce measurements of V 2 . This procedure resulted in a total 
of 35 estimations of V 2 and a summary of each night is given in 
Tabled 

It should be noted that the uncertainties in the calibrator an- 
gular diameters, when compared to the measurement uncertainty, 
have a negligible effect on the final V 2 but for consistency, were 
included in the uncertainty calculation. 



3 ANGULAR DIAMETER 

The brightness distribution of a star can be modeled, in the simplest 
case, to be a disc of uniform irradiance with angular diameter #ud • 
The theoretical response of a two aperture interferometer to such a 
model is given by 

2 

(1) 



\v\ 2 = 



2Ji(7ri6|6» UD /A) 



7r|b|6»uD/A 



where Ji is a first order Bessel function, b is the baseline projected 
onto the plane of the sky and A is the observing wavelength. For 
stars that have a compact atmosphere only small corrections are 
needed to account for monochromatic limb-darkening - these cor- 



rections can be found in lDavis. Tango &~Bo oth ( 2000 

An additional parameter, A, is included in the uniform disc 
model to account for instrumental effects arising from the differing 
spectral types of (3 Hyi and the calibrators. Equation Q} becomes 



IVf 



A 



2Ji(tt|6|6»/A) 



n\b\8/\ 



(2) 



The estimation of #ud and A was completed using \ 2 minimi- 
sation with an implementation of the Levenberg-Marquardt method 
to fit equation ([2} to all measures of V 2 . The formal uncertain- 
ties derived from the diagonal elements of the covariance matrix 
(which are calculated as part of the \ 2 minimisation) were veri- 
fied by Monte Carlo simulations - the visibility measurement er- 
rors may not strictly conform to a normal distribution and equation 
([2} is non-linear. These simulations subjected the model visibilities 
of each observation to a Monte Carlo realisation of the (assumed) 
normal measurement error distribution to produce synthetic data 
sets. Estimates of the model parameters are then found using the 
synthetic datasets, thus building a distribution of each parameter. 

The reduced \ 2 of the fit was 1.29 implying that the measure- 
ment uncertainties were underestimated. We have therefore scaled 
the uncertainties in the fit parameters by Vl-29 to obtain a reduced 
X 2 of unity. Final values of #tjd and A (with the associated la 




20 40 60 

Projected Baseline (m) 



80 



Figure 1. Nightly mean V 2 measures with the fitted uniform disc over- 
layed. Inset: all V 2 measures at the (nominal) 80 m SUSI baseline. 



parameter uncertainty) are 2.156 ± 0.017 mas and 1.008 ± 0.007 
respectively. The data, with the fitted uniform disc model overlaid, 
are shown in Figure [T] The fitted value of A is close to unity in- 
dicating that the effect of the differing spectral types of the stars 
is very small and can be neglected. Furthermore, the resolution of 
SUSI during observations was such that the calibrator stars were 
essentially unresolved. 

Analysis of the effect of the wide observing band was com- 
pleted using a G2 IV flux distribution (Davis et al. in preparation). 
Wide bandwidth effects are only significant when the interferom- 
eter's c oherent field-of-view is smaller than the angular extent of 
source jTango&Davisll2002h . The coherent field-of-view during 
observations was found to be greater than 7 mas hence bandwidth 
smearing can be considered negligible. However, the interferom- 
eter's effective wavelength when observing a G2 IV star is ap- 
proximately 696.6 ± 2.0 nm (Davis et al. in preparation). Subse- 
quently, the final fit to equation (O was completed w i th the new 
effective wavelength. Using the work of iDavis et alj (2000), the 
limb-darkening correction was determined to be 1.047 using the 
following parameters: T e « = 5872 K, logy = 3.95 and [Fe/H] = 
-0.17. These parameters (with the exception of [Fe/H]) are the re- 
sult of an iterative procedure whereby initial values, obtained from 
SIMBAD, were refined by subsequent fundamental parameter de- 
termination (see Section [4]l. We note, however, that the change in 
T e s and log g only resulted in a decrement of the last digit in the 
limb-darkening correction. The parameter [Fe/H] is the mean of 
the values listed on SIMBAD and we estimate an uncertainty in the 
limb-darkening correction of 0.002. Therefore the limb-darkened 
diameter is #ld = 2.2 57 ± 0.019 mas an d carries the caveat that 
the Kurucz models that lDavis et alj |2000) based their work upon 
are accurate for this star. 



4 DISCUSSION 

Formulae to calculate the stellar luminosity (L), effective temper- 
ature (T c ff ) and radius (R) from the observable quantities angular 
diameter (9), bolometric flux (/) and parallax (tt p ) are: 
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Table 2. Summary of observational data. The night of the observation is given in Columns 1 and 2 as a 
calendar date and a mean MJD. The nominal and mean projected baseline in units of metres is given in 
Columns 3 and 4 respectively. The mean squared visibility, standard deviation and number of observations 
during a night is given in the last three columns. 



Date 


MJD 


Nominal 
Baseline 


Mean Projected 
Baseline 


Calibrators 




V 2 


a 


#V 2 


2004 October 08 


53286.56 


60 


40.87 


e Pav, 5 Hyi, 


eHyi 


0.359 


0.019 


2 


2004 October 29 


53307.48 


5 


3.40 


e Pav, 5 Hyi 




1.037 


0.033 


4 


2005 September 20 


53633.64 


40 


27.08 


e Pav, S Hyi 




0.671 


0.021 


6 


2005 November 1 1 


53685.45 


80 


54.40 


e Pav, 8 Hyi, 


eHyi 


0.145 


0.005 


6 


2005 November 12 


53686.42 


80 


54.23 


e Pav, 5 Hyi, 


eHyi 


0.141 


0.009 


5 


2005 November 13 


53687.47 


80 


54.29 


e Pav, 8 Hyi, 


eHyi 


0.146 


0.006 


7 


2005 November 27 


53701.43 


5 


3.40 


e Pav, 5 Hyi, 


eHyi 


0.983 


0.024 


5 



^ J 1 eff 

Figure 2. Relationships between the observable quantities (inner) and fun- 
damental stellar parameters at the vertices. 
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where a is the Stefan-Boltzmann constant and C is the conversion 
from parsecs to metres. The relationship between these quantities 
can be visualised as Fig. [2] where the inner quantities are the ob- 
servables and the vertices are the fundamental stellar parameters. 

Our measurement of the angular diameter of (3 Hyi, after cor- 
recting for limb darkening, has an accuracy of 0.8%. Combining 
this with the Hipparcos parallax (0.4%) allows us to determine the 
radius of the star with an accuracy of 0.9%, as given in Table [3] 
This radius determination locates the star in the H-R diagram with 
much greater accuracy than was previously possible (see Fig. [3} 
and will be extremely important for calculating theoretical models 
for (3 Hyi. 

Our angular di ameter measurement, in combina tion with the 
bolometric flux of iBlackwell & Lvnas-Gravl dl998l) (the associ- 
ated uncertainty estimation has been taken from |Pi Benedetto! 
1998) yields an effective temperature of 5872 ± 44 K. This value, 
with an accuracy of 0.7%, is hig her than those pres ented in 
IBlackwell & Lvnas-Gravl d 19981) and lDi Benedetto! dl998l) : 5710 ± 
29 K and 5774 ± 52 K respectively. W hile all three temperature 
estima tions use the bolometric flux of IBlackwell & Lvnas-Gravl 
(1998), only our value uses a direct measurement of the angular 
diameter (and hence radius) of /3 Hyi. Furthermore, a recent cali- 
bration of the MK system for late A-, F- and early G-type stars by 
iGrav et alj d200ll) yields a value of 5850 K for a G2 IV star and the 
recent spectroscopic analysis by da S ilva et all |2006) produced an 
effective temperature of 5964 ± 70 K for (3 Hyi, both of which are 
consistent with our value. 
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3.750 



Figure 3. The 1-sigma constraints presented in this work on the effective 
temperature (vertical), luminosity (horizontal) and radius (diagonal), shown 
in an H-R diagram. 



The luminosity of (3 Hyi, found from the Hipparcos 
parallax and bolometric flux, is 3.51 ± 0.09 L© (following 
iBahcall. Pinsonneault & Basj l200ll we adopt L Q = 3.842 x 
10 26 W with an uncertainty of 0.4%). The constraints these fun- 
damental stellar values place on the location of (3 Hyi in the H-R 
diagram are shown in Fig[3] 

Calculating theoretical models for [3 Hyi is beyond the scope 
of this paper. However, we are able to combine our radius measure- 
ment with the mean density, inferred from the observed value of the 
large frequency separation in the asteroseismic data, to determine 
the mass using: 



M = -ivpR 3 . 



(6) 



Taking the value for the mean density of (3 Hyi, determined to a 
precision of 0.6% bv lBedding etall d2007h using asteroseismology 
(see Tabled, we calculate the mass to be 1.07 ± 0.03 Mq. This 
is consistent with, but much more precise than, values in the lit- 
erature that were estimated fr om modelling the p osition of [3 Hyi 
in the H-R diagram: 1 1 Mm dDravins et alj|l998l) , 1.10±°°^ M w 
dFernandes & Monteirol I2OO3I) . 1.17 ± 0.05 M dda Silva et all 
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Table 3. Physical parameters of j3 Hyi. 



Parameter 


Value 


Uncertainty (%) 


Source 


0ld (mas) 


2.257 ±0.019 


0.8 


this work 


7T P (mas) 


133.78 ±0.51 


0.4 


ESA97 


/ (fO" 9 Wm- 2 ) 


2.019 ±0.050 


2.5 


BLG98, 








D1B98 


L(L ) 


3.51 ±0.09 


2.6 


this work 


T efI (K) 


5872 ± 44 


0.7 


this work 


i?(R ) 


1.814 ±0.017 


0.9 


this work 


p(P@) 


0.1803 ±0.0011 


0.6 


BKA07 


M(M ) 


1.07 ±0.03 


2.8 


this work 


log g 


3.952 ± 0.005 


0.1 


this work 



ESA97 :lESAHl997l) ; BLG9 8: lBlackwell & Lvnas-Gra-yM 19981) : 
DiB98: lDi Benedettcl jl998l) ; BKA07: bedding et all 120071) 



2006). Finally, we can combine our radius with the mean density to 
produce an estimate of the surface gravity, 

g = ^GTipR, (7) 

where G is the universal constant of gravitation. The value we ob- 
tain for the surface gravity leads to log g — 3.952 ± 0.005, which 
has a precision of 0.1%. 



5 CONCLUSION 

We have presented the first angular diameter measurement of the 
G2 subgiant j3 Hyi. In combination with literature values for the 
bolometric flux and parallax, the angular diameter has constrained 
the stellar radius and effective temperature. The radius measure- 
ment, combined with the mean density determined from asteroseis- 
mology, allows the most accurate mass estimate of j3 Hyi to date. 
Indeed, this is perhaps the most precise mass determination of a 
solar-type star that is not in a binary system (apart from the Sun). 

The constraints on L, T e ff, R, M and \ogg that we present 
will be invaluable in the future to critically te st theoretical mod- 
els of (3 Hyi and its oscillations. As stressed by |Brown~ & Gilliland 
for example, oscillation frequencies are of most importance 
for testing evolution theories when the other fundamental stellar 
properties are well-constrained. 
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